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Abstract 

It can be fruitful to view two-component physical systems of attractive 
monomers, A and B, "chemically" in terms of a reaction A -|-B ^ C, where 
C = AB is an associated pair or complex. We show how to construct free en- 
ergies in the three-component or chemical picture which, under mass-action 
equilibration, exactly reproduce any given two-component or "physical" ther- 
modynamics. Order-by-order matching conditions and closed-form chemical 
representations reveal the freedom available to modify the A-C, B-C, and 
C-C interactions and to adjust the association constant. The theory (in the 
simpler one-component, i.e., A = B, case) is illustrated by treating a van der 
Waals fluid. 
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1. Introduction and background 

The clustering of atoms and molecules in equilibrium is a nearly ubiquitous experimen- 
tal phenomenon, occurring even among neutral, symmetric atoms like argon Not 
surprisingly, then, theoretical investigation of chemical association has remained a basic 
topic in physical chemistry since the work of Dolezalek early this century. Advances 
in microscopically-based statistical mechanical descriptions include those of Hill 0] and 
Wertheim [^. Of greater interest here, however, are the more practical "thermodynamic" 
approximations, such as those of Heidemann and Prausnitz , who proposed a closed-form 
analytical approximation embodying all cluster sizes, and of Ebeling whose approach for 
electrolytes matched exactly known virial coefficients. More recent treatments are discussed 
by Anderko [^] and Economou and Donohue [|1^ . 



The present discussion is limited to chemical formulations describing only "pairing" 
and employing a single association constant K{T) = K2{T). In the context of electrolyte 
theory — which serves to illustrate general characteristics of chemical approximations — the 
importance of dipolar pairing was first stressed by Bjerrum [|TT|, and has proved crucial 



in describing the critical region of the "restricted primitive model" [|T^|T3 • electrolytes, as 
in other systems, one expects cluster-inducing attractions to dominate at low temperatures, 
T: loosely speaking, clustering descriptions attempt to represent the partial "saturation" of 
such interactions, insofar as clustered particles are supposed to be in close contact. 

The theoretical description of the temperature variation of association presents serious 
challenges and entails basic restrictions on the form of the association constant [p!2| -[15[|. As 
T rises and excluded- volume effects dominate the thermodynamics, one expects heuristically 
the degree of clustering to drop and eventually to vanish at high T. However, with Lee 
we have recently shown that electrolyte pairing theories with association constants, K{T), 
which vanish at high T yield negative specific heats and so violate thermodynamic convexity 
requirements (resting on the Second Law of Thermodynamics |[T^) for significant regimes of 
temperature and density, p. Moreover, the underlying mechanism |T2||T^ proves not specific 



to electrolyte theories. On the other hand, acceptable forms, like the simple expression 



K{T) = K^e^°/'^ , do not vanish at high T; but the resulting finite value, K = when 
T ^ oo, then imposes some residual pairing even when repulsive forces totally dominate. 
Chemical descriptions employing appropriate expressions for K(T) are, therefore, forced to 
accommodate pairing even for hard-core particles at infinite temperature! 

A pairing description of hard-core particles, in addition to seeming rather unphysical, 
is not easy to accomplish accurately. Consider a model of attractive, spherical monomers, 
A and B, with hard-core repulsions. A fairly good "physical" equation of state for the 
monomers when T ^ oo is provided by the Carnahan-Starling (CS) theory and its extensions 
||T7| , p!8[| . However, once C = AB pairs are also considered to be present in mass-action 
equilibrium at high T, as implied by = K{T — oo) > 0, no comparably accurate 
c/iemzca/ hard-core descriptions seem available. (See [jl5| for further details.) 

These considerations lead to the basic question addressed here: "How can one incorporate 
into an associative or chemical framework as much of the exact thermodynamics as is reliably 
known, in a tractable way that does not violate thermodynamic convexity requirements?" 
In particular, one wants to maintain the accuracy of the monomeric repulsive-core physics 
without violating thermodynamic principles; but also desires to incorporate known virial 
or high-T expansions. In all cases, one hopes to benefit at low temperatures from the 
associative, chemical description which, at least at lower densities, becomes increasingly 
reahstic. (See, e.g., 

In pursuit of these goals, we have studied systems consisting of two distinct monomeric 
species, A and B, which are assumed to interact via short-ranged or van der Waals, but not 
necessarily pairwise-additive, forces. Such systems are approached in two ways: (i) via a 
"direct," "monomeric," or "physical" picture, in which the thermodynamics are given solely 
in terms of the temperature, T, and the two monomeric number densities pa and pb] and (ii) 
in terms of an "associative," "pairing," or "chemical" picture where the thermodynamics 
are specified using three densities, namely, pa and pb for the "free species," and pc for 
the associated pairs, "complexes," or "two-particle clusters," C, along with a mass action 
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law relating the equilibrium densities corresponding to the reaction A + B ^ C. Most of 
the principles developed will apply also to ionic fluids; the analysis simplifies in the one- 



component case where A = B and C = A2 ||T5 |. 



We show how to match the chemical and physical descriptions explicitly and ex- 
actly. More concretely, given a "physical" Helmholtz free energy density f(T;pa,Pb) = 
—A{T; Pa, Pb) I ksTV for a monomeric system of volume V, we demonstrate how to construct 
a chemical three-component free energy density, f+{T; pa,pb,Pc), which — when minimized 
with respect to the adjustable pair density, pc, for a given association constant, K{T) — pre- 
cisely reproduces the thermodynamics implied by /. To express this fundamental matching 
criterion analytically, let p^^ = pc(T; pa, Pb) denote the equilibrium pair density determined 
via the appropriate mass-action law [i.e., by matching chemical potentials or minimizing 
A+(T; Pa, Pb, Pc)]'- then the desired relation is 

U{T; Pa, Pb, pT) = RT; Pa + p^, Pb + P^)- (1) 

Although the concept of this precise equivalence is not new (see, e.g., P,|7|,|5|, p!P|J2Il| ] ) we are 
unaware of any similar explicit and exact results for chemical pictures employing only pairs. 

One might feel that if the physical thermodynamics are already known, seeking an equiv- 
alent pairing description is of purely academic interest and of little practical use (beyond 
surmounting the high-temperature difficulties discussed above). Our primary hope, however, 
is that the analytical prescriptions discovered may be used as tools and guides for incor- 
porating information that is known exactly or reliably (for example, in terms of some well 
understood reference system) into approximate associative thermodynamic descriptions. In- 
deed, we are applying the methods to studying the restricted primitive model of electrolytes 
and its extensions pT| . 



2. Physical and chemical thermodynamics 

The "physical" thermodynamics are given solely in terms of the monomer densities pa 
and Pb which, for clarity, we call u and v, respectively, so that the overall monomer density 
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is just p = u + V. The free energy density in d spatial dimensions can then be written 

/(T; v) = u[l - In (Af n)] + - In (Afi;)] - F(T; n, v), (2) 

where we suppose that the excess free-energy density, F, embodies the "true thermodynam- 
ics" which the chemical picture must match. As usual, Aa(T) and Af,(T) denote the standard 
thermal de Broglie wavelengths. For the most part, the temperature dependences will not 
be displayed since the basic manipulations will be performed at fixed T. 

In the "chemical" or associative picture, the overall monomer density is now 
p = u + V + 2w, where w = p^ is the pair or cluster density. Correspondingly, we intro- 
duce the augmented or chemical (reduced) Helmholtz free-energy density via 

/+(T; u,v,w) = u[l-\n (Afw)] +v[l-\ii (A^t;)] +w{l-\n [Af w;/ir(T)]} 

+ fl-{T;u,v,w), (3) 

where, to ensure consistency with the classical momentum integrals in the direct partition 



function, we suppose A^ = y/A^JVb [|14|. Notice that the association constant, K{T), appears 



here in the role of the internal partition function of an AB = C pair [|T^|T^^; for the 
present purposes, however, K{T) may be supposed arbitrary. 

Now, the excess chemical free-energy density may be decomposed according to 

/^^(■u, f , w) = —F{u, v) — uvE{u, v) — wD(u, v, w), (4) 

where the first term, —F{u, v), represents just the physical-picture excess free energy which 
must be recovered in the absence of association, i.e., when w = pc —>■ 0. Then, the "aug- 
mented monomer-interaction function," E{u,v), represents the modifications to the A-B 
interactions inherent in the association process Q]. Finally, the "cluster- interaction func- 
tion," D{u,v,w), embodies the cluster-monomer, A-C and B-C, and cluster-cluster, C-C, 
interactions that must be introduced. 

The condition of chemical equilibrium under the reaction A + B ^ C follows in standard 
fashion [by equating the chemical potentials or by minimizing — /_|_(T; u, v, w)] leading to the 
mass-action law, 



w = Kuv-ia-ihhc = Kuve^^''''''''\ (5) 

where the activity coefficients, 7o-, and the excess chemical potential difference, M{u,v,w), 
are given by 

^-^-^j/TK^'^)- (6) 
Of course, (j^) represents a highly nonlinear equation which is normally quite intractable. 



3. Matching via virial expansions 

The free energies for both the physical and chemical pictures may be expanded in powers 
of monomer and pair densities which permits the establishment of order-by-order matching 



conditions. (Compare with |19|.) One may iteratively solve the mass action relation 
to obtain the solution w = Wc^iu, v) formally to any desired order. An expansion of the 
fundamental matching relation (0) solely in terms of u and v may then be obtained by the 
substitution w =^ Weq. The nth order matching conditions result from equating coefficients 
of u'^v"'^'' (for all permissible q) in the expansion ([l|). 

Letting Fim, Eim, and Dimk represent the coefficients of u'f™ in F and uvE, and of 
^lyrn^k ^jj^ respectively [see (|)], we find, for the two lowest orders, 

n = 2 : En = K, (7) 

n = 3: E2i + KDioi = lK'' + Ki2F2o + Fu), (8) 

Eu + KDon = + K{2Fo2 + Fii). (9) 

The higher-order relations for both the two- and single- component (A = B) cases have a 



completely analogous structure and have been obtained explicitly up to n = 5 ||T5|. 

These relations reveal physically significant features of the chemical picture. The lowest 
order relation (|^) — which is implicit, if not more-or-less explicit in much previous work. 



e.g. [P,[7|, p!l| , p!U| , P^p5| — admits no adjustable parameters, assuming K is specified. Since K 
must be positive, it implies that the chemical A-B second virial coefficient (the coefficient 



of uv in p/ksT) Buo = Fu + En, increases relative to its "bare" value, Bu = Fu, in the 
physical picture. In other words, the original A-B interactions must become more repulsive 
in the chemical picture. 

In next order the relations reveal the freedom implicit in the formulation: thus in (^, 
only the combination {E21 + KD iqi) is constrained to match the physically specified, lower- 
order terms on the right-hand side. Consequently, the A-C monomer-cluster interaction, 
which is "encoded" in the second virial coefficient -Dioi, can be chosen arbitrarily; likewise, 
the B-C interaction. In nth order on finds similarly, that one may choose all the 

cluster interaction coefficients Dimk (and hence the C-C interactions, the A-B-C three-body 
coupling, etc.) arbitrarily by fixing the corresponding modifications of the bare monomer 
interactions as embodied in the n — 1 coefficients Eim with l,m > 1, I + m = n [p!5[] . 



4. Exact closed-form chemical formulation 

The latitude discovered in the chemical formulations may, it transpires, be exploited to 
construct exact, closed- form chemical representations for arbitrary physical thermodynam- 
ics. In particular, the solution of the intractable mass-action relation may be side-stepped 
at a comparatively small cost in fiexibility |]TB[. To this end, let us focus on the excess chem- 
ical potential difference, M{u,v,w), given by (^). If the mass-action equation (|^) had been 
solved for w = Weq{u,v) — given the interaction functions F, E, and D in (|) — one could 
explicitly find M as a function only of u and v. This suggests that we may choose, at our 
discretion, a "solved" chemical potential difference function Mo{u,v) [which then represents 
the exact equilbrium value of M{u, v, w)] and impose the constraint 

{du + dv) [F{u, v) + uvE{u, v) + wD{u, v, w)]^^ - [wD{u, v, w)]^^ = Mo{u, v). (10) 

Here we adopt = d/dx and use the subscript "eq" to indicate that one must set w = 
Weq{u, v) after differentiation. But Weq may now be obtained trivially: by comparing (|10|) to 
(|) and (^), one sees that the mass-action law has, indeed, been solved in closed form giving 
just 
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Weq{u,v) = Kuv exp[Mo{u,v)]. (11) 

Now the implications of the constraint (0) are not so obvious; however, it is hnear in 
Mq, F, E, and D and, consequently, one finds |[T5| that much freedom remains explicitly 
available. In particular, it proves convenient to specify the cluster-interaction function, D, 
by choosing the "component functions," Di and D2, in the form 

D{u, V, w) = Do{u, v)Di{u, V, w) + D2{u, v, w). (12) 

Physically reasonable and effective choices for Di and D2 are discussed below. The remaining 
"coefficient" Dq{u,v), and the augmented monomer interaction function, E{u,v), can then 
be expressed simply in terms of F, Mq, Di, and D2: see the Appendix. 

Once Dq{u,v) and E{u,v) are calculated via the straightforward prescription in the 
Appendix, the exact chemical description is complete: minimizing the free energy (^ with 
respect to w (i.e., pc) should now precisely reproduce the original thermodynamics embodied 
in F{u,v), the physical free energy. 

5. Faithfulness of the exact representations 

Although the closed-form solutions for the chemical free energies constructed in the 
previous section are formally exact, one can ask: "Do they actually work?" To answer 
this question, consider, for simplicity, a one-component system (A = B) undergoing the 
reaction 2 A ^ C (although similar considerations apply to the A + B ^ C system). Given 
a chemical free energy, f^{u,w), the task of finding equilibrium solutions Meq(p) and Weq(p) 
for a fixed overall monomer density p = u + 2w, can be phrased as seeking the minimum 
chemical free energy at fixed p. Now the exact chemical formulation in Sec. 4 guarantees 
that the mass-action law (^ is satisfied (with v = u); but that ensures only that f^{u,w) 
has an extremum on the required locus w = Ku'^e^'^^^^^ [see (pT])]. It is not guaranteed, 
however, that this extremum is a free-energy minimum, let alone the global minimum, as 
essential for a sensible chemical description. 
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Thus one may find that not every representation is "faithful" — i.e., exactly reproduces 
the original, physical thermodynamics — for all temperatures and densities. One may expect 
that the desired physical solution will be the global free-energy minimum (for arbitrary Di 
and D2, say) when the densities p, w, u (and v) are sufficiently small; but failures can, 
indeed, arise at higher densities for inadequate choices of D. 

To illustrate this, consider a one-dimensional fluid of hard rods of length h for which the 
physical "excluded-volume" excess free energy is 

F^^{u) = -u\Yi{l-hu). (13) 



(See, e.g., [QO To obtain an exact chemical picture with the choice of association constant, 
say, K = b, one may try first the "minimal representation" generated by the simplest 
choice, namely, Mq = 0, -Di = 1, and = in (0) and (0), with Dq{u) then determined 
by the one-component versions of (Al)-(A7) Numerical and analytical investigations 
reveal, however, that this "minimal" chemical picture remains faithful only up to p* = pb = 
p/Pmax — 0.26 for larger densities the global chemical free-energy minimum becomes 
totally unphysical! 

Nevertheless, we have found that the situation can be easily remedied by using physical 
intuition to select a cluster interaction function Di that refiects the hard-core aspects of the 
A-C and C-C interactions that one would reasonably expect. Thus, in parallel to ([T3|) , the 
"excluded-volume" specification 

HC : wDi = —{u + w) In {1 — buU — b^w) + u\n{l — buu), (14) 

along with D2 = Mq = (as previously) proves very effective. Explicitly, the natural choice 
bu = b, bw = 2b is found, by numerical examination, to yield faithful representations up to 
p* > 0.98 for assignments of association constants, K, ranging from 10~^b up to 10^6. 

Furthermore, fully comparable results are obtained for {d = l)-dimensional hard-rod 



A+B mixtures and, in = 3 dimensions, for the corresponding Carnahan-Starling p7|JT8 
hard-sphere "physical" thermodynamics (regarded as exact) |T^. To test and illustrate the 
theory further we consider systems with attractive interactions. 
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6. Application to a van der Waals fluid 

The van der Waals (vdW) equation gives a surprisingly good semiquantitative description 
of gas-liquid criticality in simple one-component real fluids which, furthermore, is exact in 



an appropriate inflnite-range limit |Q. The corresponding excess free energy is 



F(T,p) = -pln(l - hp) - h{e/kBT)p\ (15) 

with 6, e > so that convenient reduced units are p* = bp and T* = kBT/e. The Boyle 
temperature, Tb, and critical parameters are then = 1, T* = ^, p* = |. See Fig. 1 for 
the coexistence curve. 

The assignment of an appropriate association constant now calls for some study m 



The choice K'^'^'^lT) = 6exp (1/T*) matches the temperature dependence of the vdW sec- 
ond virial coefficient at high T — compare with (|^) and Refs. [[F|,P,p!2|- [l^j2^ — and also 



satisfies convexity criteria needed (in approximate representations) to avoid the Second-Law 
violations referred to above |T2| , p!5| . It has been adopted for the tests reported here. 

Chemical representations of the vdW thermodynamics (|1^) can be assessed in terms of 
faithfulness boundaries, or loci: for every temperature, T, there will be a density, say po{T), 
up to which a given representation is faithful. As reported for the hard-core models, the 
minimal formulation (Mq = D2 = 0, Di = 1) fares poorly: see plot (a) in Fig. 1. The 
results are strongly temperature-dependent below T* ^ 0.4 and the critical point lies well 
outside the faithful domain! If one uses the temperature-independent "excluded-volume" 
cluster-interaction function Di that that accounts for A-C and C-C repulsions and proved 
so successful for hard rods — namely, (0) with b^ = ^b^^ = b — one certainly improves the 
behavior at high T*: see Fig. 1, plot (b); but below little is gained. 

To obtain better results for low T*, it is reasonable to introduce attractive (i.e., negative) 
terms into D. Embodying these for convenience only in D2{u,v,w), we found that the 
natural form in which K[T) regulates the strength of the attractions, namely, 

HCK : wD2{u, w) = -{d2iuw + d22w'^)K{T). (16) 
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proves successful. As evident from plot (c) in Fig. 1, imposing only A-C or monomer- 
cluster interactions, via the choice = 2, ^22 = 0, betters the pure HC representation 
significantly. Increasing the A-C coefficient to ^21 = 3 yields faithfulness for the full set 
of tested temperatures in the range 0.05 < T* < 00 and densities 0.01 < p* < 0.99. The 
simple choice ^21 = ^22 = 1 with both A-C and C-C attractions also proves faithful for 
the entire test set. This last formulation is probably a reasonable starting point for the 
chemical description of any model with a choice of association constant that does justice to 



the monomer attractions at low T 15 



7. Concluding remarks 

The exact chemical representations presented above permit — for the first time — pre- 
cise thermodynamic characterization of the interactions which must arise in an associative, 
A -|- B ^ C description of a system of A and B monomers. While monomer-monomer inter- 
actions are modified repulsively (since "paired" configurations are not permitted ) success- 
ful, i.e., faithful, chemical prescriptions explicitly incorporate both repulsive and attractive 
monomer-pair and pair-pair interactions. The use of just one representative or dominant 
type of cluster (i.e., pairs) clearly affects the nature of the required interactions: the pairs 
must, in effect, stand in for clusters of all sizes and, at higher densities, may therefore be 
viewed as "renormalized," "dressed," or "solvated" as discussed in ||15|. Except at low den- 
sities, therefore, pairs cannot act simply as literal two-body objects if the true "physical" or 
monomeric thermodynamics are to be accurately described. 

We hope that the exact chemical formulations will prove useful in developing improved 
associative approximations for systems with attractive forces. The explicit chemical repre- 
sentations found for hard-core reference systems should be of direct use in combination with 
existing associative approximations for attractions. Enforcing the matching relations of Sec. 
3, moreover, might further refine such theories. On a slightly different tack — starting from 
scratch and drawing solely on the "thermodynamic" conception of pairs presented here — 
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a worthwhile minimal approximation scheme might combine (i) an exact chemical formu- 
lation for a reference thermodynamics with (ii) known, low-order virial coefficients as the 
only further input. But regardless of the success of particular formulations, knowledge of 
the exact chemical expressions for simple attractive systems — which reproduce the known 
thermodynamics of, say, one-dimensional models or a van der Waals fluid — may usefully 
inform the construction of approximate chemical free energies for less tractable systems. 
Applications to electrolyte models are underway |2l[ . 
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APPENDIX A: EXPLICIT CHEMICAL REPRESENTATIONS 



The following expressions complete the chemical representation embodied in (|T^ 
and (IT 



Do{u,v) = [L{u,v) - d^{wD2)]/ [d^{wDi)] 



icq 



L[u,v\Mo,F] = [j{u,v) + {u + v)K{u,v)y Q{u,v), 

Q{u, v) = l + Kuve''^ + ^] + K{u + t;)e^° 



J(m, v) 
K{u, v) 



du dv 



OF OF fdK diC 

du dv \ du dv 



Mo{u,v) 



K[l - Mo{u, t;)]e*^° + S(u, v; fsTe*^") + AF(u, v; Kuve^^°) 
S{u, v; x) = u~^{l + ux) In (1 + ux) + t>~"'^(l + vx) In (1 + vx) — 2x, 
AF{u, v; x) = [F{u + uvx, v + uvx) — F{u, v)]/uv. 



(Al) 
(A2) 
(A3) 

(A4) 
(A5) 
(A6) 
(A7) 



The analogous results for the 2A ^ C reaction are given in Ref. |TH . 
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FIGURES 
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FIG. 1. Density-temperature plane for a van der Waals fluid with loci indicating the faithfulness 
boundaries for various chemical representations: (a) minimal, (b) "HC," and (c) "HCK"; see text 
for details. Note the nonlinear temperature scale above the Boyle temperature at = 1. 
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